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Abstract 

We evaluate theoretical uncertainties to the method of Pais and Treiman to 
measure the tttt phase shifts in K e t± decays. We find that they are very small, 
below 1%. 
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1 Introduction 



The measurement of the tttt phase shifts near threshold is a key issue of low en- 
ergy hadronic physics. First, these phases enter the phenomeno logical analysis of 
many different low energy scattering processes, weak decays and, in particular, CP- 
violating K-decays. A model independent determination of irir phases would consid- 
erably improve our understanding of the corresponding hadronic matrix elements. 
Furthermore, Chiral Perturbation Theory (CHPT) the systematic low energy 
expansion of QCD amplitudes, provides a link between the low energy mr scatter- 
ing data and the non perturbative chiral structure of the QCD ground state: The 
standard CHPT has a clean prediction for these quantities [TJ, [|, and would not 
be able to explain a large discrepancy with experimental data. On the other hand, 
the proposed generalization M of the usual low energy expansion, called generalized 
CHPT, expects a somewhat stronger 1=0 S-wave nir interaction. Within the latter 
scheme, a discrepancy between the standard CHPT and experiment would be inter- 
preted as a manifestation of unusually low values of the quark-antiquark condensate 
(qq) and of the ratio of strange to non strange current quark masses. 

It is known since a long time that the main source of model independent ex- 
perimental information on low energy nir phases comes from K ei decays. In these 
decays, a complete im phase-shift analysis could be performed in principle, assuming 
nothing more than unitarity or the Watson final state interaction theorem. How- 
ever, such a complete procedure would require a detailed amplitude analysis of K e4: 
to be performed with respect to all five kinematical variables. To avoid this task, 
which is rather problematic in practice, Pais and Treiman have proposed || to mea- 
sure 7T7T phases in a much simpler way which, in addition to the Watson theorem, 
assumes that higher partial waves beyond the P-wave are small. A variant of the 
Pais- Treiman method has already been used (together with additional assumptions) 
in a 1977 high statistics experiment by Rosselet et al. @, but the data obtained 
still show very large error bars. Newly planned experiments like KLOE at DA$NE, 
will hopefully be able to reduce the errors sizeably enough to decide between the 
theoretical alternatives mentioned above. In view of this improvement on the exper- 
imental side, it is worth to check what kind of uncertainties affect the Pais- Treiman 
method from the theoretical point of view. This is the purpose of this note. 

2 Kinematics and form factors 

We discuss the decay 

K + (p) - n + ( Pl ) 7T-(p 2 ) e + ( Pl ) u e ( Pu ) , (2.1) 

(2.2) 

and its charge conjugate mode. We do not consider isospin violating contributions 
and correspondingly set m u = rrid = rh, a QED = 0. 
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The full kinematics of this decay requires five variables. We will use the ones 
introduced by Cabibbo and Maksymowicz 0. It is convenient to consider three 
reference frames, namely the K + rest system (Ex), the 7r + 7r~ center-of-mass system 
(E27r) and the e + v e center-of-mass system Then the variables are s^, the 

effective mass squared of the dipion system, si, the effective mass squared of the 
dilepton system, 9 n , the angle of the n + in E 27r with respect to the dipion line of 
flight in Ex, 9i, the angle of the l + in E^ with respect to the dilepton line of flight 
in Ex, and 0, the angle between the plane formed by the pions in Ex and the 
corresponding plane formed by the dileptons. The range of the variables is 



2 < (Mx - m,) 2 

2 ^ _ /_ , „ \2 ^ (■*, /—\2 



iM 2 < s n = (p! + p 2 ) 2 < (Mx - m z 



mf < Sl = (pi + p v y < (Mx - v^; 

< 9 n ,9i <tt,0 < <p < 2?r. (2.3) 

It is useful to furthermore introduce the following combinations of four vectors 

P = Pi+p 2 , Q = Pi~P2, L = pi+p u . (2.4) 

Below we will also use the variables 

t= (p 1 -p) 2 ,u= (p 2 -p) 2 ,v = t-u. (2.5) 

These are related to s n , si and 6 n by 

t + u = 2M 2 + M 2 K + si-s n , 

v = -2a v XcoB6 v , (2.6) 



where 



<r w = (1-4M 2 / S?r )2 , 
X = -A^M 2 ,^) 



X(x,y,z) = x 2 + y 2 + z 2 — 2(xy + xz + yz) . (2.7) 
The matrix element is (we always neglect the mass of the electron in what follows) 

T = ^Lv: s u(Pu)l,il - IsHPiW ~ ^) , (2.8) 



where 



h = < 7r+( Pl )7r-(p 2 )out I 1^(0) I K + (p) >; I = V,A , 
A, = -t-^[P,F + Q,G + L,R] , (2.9) 

Mx 
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and e i23 = 1- 

The form factors F, G, R and H are analytic functions of the variables s n , t and 
u. The partial decay rate for (|2.1|) can be expressed as 



dTc, — G 2 p I V, L 1 



Vus I 2 N(s n ,si)J 5 (s n , s h 9 n , 8 h (j))ds n dsid(cos 9 7T )d(cos9i)d<p , 
N(s^ Sl ) = a 7T X/(2 13 ir Q M 5 K ) , (2.10) 

where 

J 5 = 2 l 

+ J 5 sin 9i ■ cos 4> + Iq 

+ Jg sin 2 9i ■ sin ! 



with 



and 



l\ + I 2 cos 29[ + I 3 sin 2 9\ ■ cos 20 + I4 sin 29i ■ cos 

• ™ c ^ -J- ^ cos 9i + ^7 sin 9\ ■ sin + Jg sin 20; • sin ■ 



h = 
h = 
h = 



;{l^i| 2 + 5(l^| 2 + |F, 
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3 1 J sin 6* 

+ 1^3 



sin 2 On 



{|F 2 | 2 -|F 3 | 2 }sin 2 ^ , 

Re(F*F 2 )sin^ 

Re^i^sin^ 
Re(F 2 *F 3 ) sin 2 9, 
Im(F*F 2 )sin^ 



y 7T ; 
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h = ~ 
2 

/s = - 

h = ~ 

I? = ~ 

h = \ Im(F*F 3 )sin^ , 

J 9 = -- Im(F 2 *F 3 ) sin 2 9 n , 



'7T 5 



(2.11) 



F 3 



X-F + a^(PL)cos^-G' 

) 1/2 G , 



Ml ' 



(2.12) 



The definition of Fx, F 2 and _F 3 corresponds to the combinations used by Pais and 
Treiman || . The form factors Ix , ■ ■ ■ , I9 agree with the expressions given in || if 
one neglects the mass of the electron. 

The form factors may be written in a partial wave expansion in the variable 9 n . 
Suppressing isospin indices, one has for F and G M 



1=0 



iSi 



cosfl^G , 
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oo 

G = J2 P !( cosd «)9ie i51 , 
1=1 



(2.13) 



where 

Pl(z) = j-Pii*) ■ (2-14) 

The partial wave amplitudes /; and gi depend on s n and si and are real in the 
physical region of K i4 decay (in our overall phase convention). The phases Si coincide 
with the phase shifts in elastic kit scattering. Accordingly, the form factors Fi and 
F 2 have the expansions: 

oo 

F 1 = X^P^^)/^ , 

1=0 

oo 

-sin^F 2 = ^(v^/^P/Vse^/' , (2.15) 

i=i 

where Pj^ are the associated Legendre functions. 



3 Evaluation of the corrections 

The method suggested by Pais and Treiman to measure the mr phase shifts in K ei 
decays is very simple and clean. It is based on the observation that the dependence of 
the differential decay rate on two of the five variables can be worked out analytically 
under very few general assumptions, as we have seen in the previous section. 

In case one neglects all the waves higher than S and P, I4 and I7 show a very 
simple dependence on the phases So and Si. [The same dependence appears in J 5 
and Is which, however, contain the anomalous form factor H, and are kinematically 
much more suppressed.] Their integral over cos^ is: 

/l 7T 
d cos 0,^4 = -Xa v {s v si) 1/2 hgi cos(<5 - S ± ) , 

cos 6J 7 = -Xo- 7r (s 7r s i ) 1/2 /ofi'isin(5o - Si) . (3.1) 

By measuring the ratio I7/2I4 one has then direct access to tan(5 — Si). To exper- 
imentally select I7 and I4 one could use different methods (fit the distribution in 
cos#; and <j), measure asymmetries, measure moments . . . ). Choosing for example 
to measure the appropriate moments, one would have, all in all, to integrate the 
distribution over four of the five variables with the weights sin0, then cos cos 9i, 
and then take the ratio. The measurement of tan(5 — Si) would be as simple as 
that. 

While on the experimental side there might be detector-dependent problems that 
could force experimentalists to adapt the Pais- Treiman method to their own peculiar 
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situation (see for example Ref. ||), it is clear from the theoretical point of view 
that the beauty of the method is based on neglecting the higher waves. It is the 
main purpose of this note to try to estimate corrections coming from higher waves 
to the Pais-Treiman formula: 

i = tan(5 - 5 t ) . (3.2) 

For that purpose we shall use the form factors predicted by CHPT at the one loop 
level. 

We start by applying the chiral power counting to the partial waves of the form 
factors. In general, for any partial wave we may write an expansion in the following 

way: 

X, = + X? > + Af > + ...}; X = /, g , (3.3) 

where the upper index stands for powers of energies or meson masses. From pub- 



lished calculations of F and G to one- loop B, 10 1 we may easily get that only three 



partial waves start at order 0(E°), while all the others start at order 0(E 2 ): 

Mu 



fo 



^{l + ff + 0(F*)} 



fi = ^{// 2) +0(^)};/>2 



9k 



-^-{8» + 0(B*)} ;*>2 . (3.4) 



Using this chiral power counting is very easy to give a "corrected" Pais-Treiman 
formula, which is accurate up to and including order 0(E 2 ): 

±L= t{m (S -6 1 ){l+AW + O(E 4 )} , (3.5) 

A (2) ( sin^o COS £p \ ~ (2) 

A ~ T ~Tx T\ 7x T\ 2^ A 02k+ig 2 k+i 

\sin(do - di) cos(d - di)/ ^ 

/ sin S 1 cosSx \ a n PL ~ (2) 

+ [sin(6 -6 1 ) oos(do -^J X jL> Al *** 

~ f^ + ^lE^n/Jf , (3.6) 
\sin(do-di) cos(d -di)/ ^ 
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where 

A lk = - f 1 dcos^P z (cos^)P fc (1) (cos^) , (3.7) 

7T J -I 

and where all the phases of waves higher than the P have been put to zero (this 
is again consistent at the order at which we are working). The ratios of sines and 
cosines start at order 0(E°), and will have also contributions of higher order, that 
we neglect at this level of accuracy. We may thus use the leading order CHPT 
expressions: 

60 = ^^^{ 2 ^- M " + ^ M ")+°( EA ) (3-8) 
*i = ^rf^ + O^) ■ ( 3 - 9 ) 

The standard CHPT predicts e = 0, whereas in generalized CHPT e is an arbitrary 
parameter, < e < 1, related to the quark mass ratio r = m s /m Eq. (3.6) then 
becomes 



fc=i 



'As n + 5M2 + lheMl\ a^PL ^ A (2) 



% 5s, + Ml + 15eM2 / X 
'6 S7r -3M2 + 15eM2 



5s n + Ml + 15eM2 



fc=i 

(2) 

i=i 



Before evaluating numerically the correction we would like to stress that the low 
energy expansion has been performed only inside the braces of eq. (|3.5|) . 

The explicit calculation shows that both in standard and in generalized CHPT 
AW is numerically very small over the whole phase space. It barely reaches 0.5 %. 
The dependence on si is rather weak, and there is no difficulty in repeating the same 
analysis for the ratio of I7 and 2/4 not treated as functions of si but averaged over 
it (as suggested by Pais and Treiman). At this order in the low energy expansion 
the correction turns out to be just (AW), where with the symbol () we mean the 
following average: 



(A) = dstX^Afr) (3.11) 



'1 



The plot of (AW) function of is shown in fig. 1. 

Since the corrections are so small, we would like to understand whether this 
is just a typical size of higher waves, or an effect depending on some miraculous 
cancellation. 
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Since F and G are dimensionless functions of s„, si and u, and depend on cos6 w 
only through is, their higher waves projections are suppressed by powers of the 
kinematical factor: 

(3.12) 



a. 



X 



Ml 



'Ml 



< 0.15 



where p w and pi stand for the momenta of individual pions and of the dilepton, 
respectively, in the dipion center of mass frame. The bound in (|3.12|) holds over the 
whole phase space. 

The dominant contribution to the sum A 1 - 2 ), eq. ( |3.10 ), comes from the D- waves 



('2') i ' i 

g\ and f^' ■ The latter contains contributions both from the F and G form factors: 



A2) 



/ 2 (2) = fT> + 2 



(2) 



a. 



PL 



X 



-92 



(2) 



(3.13) 



where f^ stands for the D-wave of F alone. Making explicit the factors ( |3.12 ), we 
may define: 

2 



(2) 



(2) _ 
92 = 



a 



X 



Ml 



ft 



(2) 



rX. 



Ml 



92 



(2) 



'(2) 



where fd ' and are expected to be smooth functions of s n and of sj. 

D- waves are a loop effect, we may guess that fd and g^ 2 -* will be of the order 
Ml / (I6n 2 F%) . An explicit calculation shows that this is indeed the case and, for a 
possible later use, we may even write the following simplified expressions for them: 
Standard CHPT gives 



(3.14) 



Since 



(2) 



92 



(2) 



Ml | 
16tt 2 F 2 

Ml 
16tt 2 F 2 ' 



-0.28) 
-0.62) 



1 + 0.12 



1 + 0.08 




Q 



(3.15) 



where 



g 2 = ( S7r -4M 2 )/4M 2 . (3.16) 

~(2) 

The above evaluation of fd is parameter free. On the other hand, a substancial 
part of the contribution to g~2^ comes from the low energy constant L 3 + 4L 2 , whose 
determination brings in some uncertainty [11|: The corresponding error in (3.15) is 
however smaller than 50%. Finally, A^ 2 ^ can be expressed in terms of the smooth 

functions fd ^ and g~2^'- 




(a) 3 / 2s* - Ml + 5eM 2 \ 
8 \5 S7T + M* + l5eM* J 
3/ g , + 3M 2 + 56M 2 W 2 (PL) ^ (2) 
2 \5s w + MZ + 15eM%) Ml 92 + "' {C ' ] 
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where the ellipsis stands for / > 2 waves. This expression gives clearly account of the 
size of the corrections. The second term is dominant with respect to the first one, 
and its typical scale is ex 2 (PL)/ (167r 2 F 2 ), which, after averaging over sj, is at most 
1.6% over the whole phase space. The remaining coefficients reduce this number by 

a factor 4. Moreover, we have calculated the dependence of fd and g^ 2 ^ on m s /rh, 

— '(2) 

in generalized CHPT[]. We have found that fd is affected by no more than 10%, 
whereas the form factor g^ 2 -* may be modified by at most a factor one half, further 
reducing (A*- 2 )). 

Higher orders will certainly modify the above numbers, but not the orders of 
magnitude, which are essentially given by the kinematics. In a sense, the limited 
effect of the uncertainties in the values of the low energy constants Lj may be 
considered as an indication in this direction. Furthermore, since the generalized 
CHPT differs from the standard CHPT by a different ordering in the perturbative 
expansion (that is, at every finite order it includes terms which in the standard 

scheme are relegated to higher orders), the weak dependence of fd and of g^ 2 ** 
on m s /m is a further indication that it is very unlikely that higher orders would 
overwhelm the strong kinematical suppression of higher waves. 

Our conclusion is then that the Pais-Treiman formula (|3.2|) is free of corrections 
up to the percent level, over the whole accessible range of energy. 
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Figure Caption 



Fig. 1: (A^ 2 )), as a function of s n . The plot refers to the standard CHPT case: 
it is calculated using eq. Q3.1C ) with e = and the F and G form factors given by 
@, [KJ, with the central values L 2 = 1.35 x 10" 3 , L 3 = -3.5 x 10~ 3 of Jll|. The 
same curve can be obtained using eq. (|3.17| ) with e = 0, and (|3.15|) . 
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